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Particle decay in the early universe: predictions for 21 cm 
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ABSTRACT 

The influence of ultra-high energy cosmic rays (UHECRs) and decaying dark matter 
particles on the emission and absorption characteristics of neutral hydrogen in 21 cm 
at redshifts z = 10 — 50 is considered. In presence of UHECRs 21 cm can be seen in ab- 
sorption with the brightness temperature Tb = — (5 — 10) mK in the range z = 10 — 30. 
Decayng particles can stimulate a 21 cm signal in emission with Tf, ~ 50 — 60 mK at 
z = 50, and Ti, ~ 10 mK at z ~ 20. Characteristics of the fluctuations of the bright- 
ness temperature, in particular, its power spectrum are also calculated. The maps of 
the power spectrum of the brightness temperature on the plane wavenumber-redshift 
are shown to be sensitive to the parameters of UHECRs and decaying dark matter. 
Observational possibilities to detect manifestations of UHECRs and/or decaying par- 
ticles in 21 cm with the future radio telescopes (LOFAR, 21CMA and SKA), and to 
distinguish contributions from them are briefly discussed. 

Key words: early Universe - cosmology:theory - dark matter - diffuse radiation. 



1 INTRODUCTION 



At 2 -^ 1000 the Universe enters the "dark ages" epoch, 
the electrons and protons recombine, and gas remains neu- 
tral until the first luminous objects emerge at 2 ~ 30 — 20. 
Neutral gas can be observed in a redshifted 21 cm line of 
neutral hydrogen in emission or absorption against the cos- 
mic microwave background (CMB). This gives a possibility 
for studying the processes associated with transition of the 
neutral universe into a fully ionized state, and for identifi- 
cation of the sources responsible for reionization (Hogan & 
Rees 1979, Madau et al 1997). At present observations of 21 
cm line are anticipated as a promising tool for diagnostics 
of the universe in the end of "dark ages" (Madau et al 1997, 
Tozzi et al 2000, Ciardi & Madau 2003). 

Possible sources of photons which can be important for 
reionization form two different groups. The first, connected 
with conventional baryon physics, involves the baryons pro- 
cessed in stellar nucleosynthesis (Shapiro & Giroux 1987, 
Miralda-Escude & Rees 1994, Tegmark et al 1994, Cen 2003, 
Ciardi et al 2003, Choudhury & Ferrara 2006) and in shock 
waves near black holes (Madau et al 1999, Oh 2001, Ri- 
cotti & Ostriker 2004). The other is connected with unsta- 
ble dark matter and can include massive neutrinos (Sciama 
1982, Scott et al 1991), superheavy X-particles (Berezin- 
sky et al 1997, Kuzmin & Rubakov, 1998, Birkel & Sarkar, 
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1998). The first luminous objects are commonly thought to 
be the principal source of the reionization. They heat gas 
in the universe through ionization by ultraviolet (UV) and 
X-ray photons. This inevitably affects the emissivity of gas 
in the 21 cm line, because the hydrogen spin temperature 
Ts depends on the gas kinetic temperature Tk, and thus ob- 
served intensity imprints the effects from the objects of the 
first generation (Madau et al 1997, Tozzi et al 2000, Ciardi 
& Madau 2003, Loeb & Zaldarriaga 2004, Zaldarriaga et al 
2004, Chen & Miralda-Escude 2004, Sethi 2005). 

Production of copious number of ionizing photons dur- 
ing the dark ages can be connected also with the origin of 
ultra-high energy cosmic ray (UHECRs) if they form from 
decaying superheavy dark matter (SHDM) particles with 
masses Mx ~ lO^'^ GeV in the so-called top-to-bottom sce- 
nario (Berezinsky et al 1997, Kuzmin & Rubakov, 1998, 
Birkel & Sarkar, 1998). The associated production of UV 
photons can have strong infiuence on cosmological recombi- 
nation (Doroshkevich & Naselsky 2002, Doroshkevich et al 
2003). 

Decaying dark matter particles, such as massive neutri- 
nos, can also contribute significantly to reionization (Sciama 
1982, Dodelson & Jubas 1994). Initial polarisation measure- 
ments of the CMB by the Wilkinson Microwave Anisotropy 
Probe (WMAP) sattelite (Spergel et al 2003) imparted a 
new impulse to this possibility. The obtained relatively large 
optical depth of the universe r ~ 0.16 suggested unrealis- 
tically strong constraints on properties of the first stellar 
objects (e.g. Cen 2003, Wyithe & Loeb 2003, see also an 
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alternative discussion in Tumlinson et al 2004), which lead 
some to consider that decaying dark matter can be at least 
a complementary source of reionization (Hansen & Haiman 
2004, Chen & Kamionkowski 2004 (hereafter CK), Kasuya 
et al 2004, Kasuya & Kawasaki 2004, Pierpaoh 2004, Mapelh 
et al 2006). The corresponding heating can change charac- 
teristics of 21 cm in emission and absorption. Very recently, 
the influence of dark matter decay and annihilation on the 21 
cm line from dark ages has been also considered by Furlan- 
etto et al. (2006a) for long living particles. For the models 
with similar parameters of decaying particles they reached 
conclusions close to ours. 

In this paper we study the effects of UHECRs and 
decaying particles on cosmological 21 cm background. 
We assume a ACDM cosmology with the parameters 
(fio, f^A, f^m, 0.b, h) = (1.0, 0.76, 0.24, 0.041, 0.73) (Spergel et 
al 2006). 



2 SPIN TEMPERATURE 

The two processes: atomic collisions and scattering of UV 
photons, couple the HI spin temperature and the gas kinetic 
temperature (Wouthuysen 1952, Field 1958) 



r. 



Tcmb + VaTk + TJcTk 



(1) 



here Tcmb is the CMB temperature, yd/a are the functions 
determined by the coUisional excitations and the intensity 
of the UV resonant photons 
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AioTk 



yc = 



CwT, 
A-ioTk 



(2) 



r, = 0.0682 K is the hyperfine energy splitting, Aio = 
2.87 X 10~^^ s~^ is the spontaneous emission rate of the 
hyperfine transition, Cio = kioun + 7en,e is the coUisional 
de-excitation rate by hydrogen atoms and electrons, the rate 
by protons is negligible, for fcio we use the approximation by 
Kuhlen et al (2006), for 7e we take the approximation from 
Liszt (2001), Pio is the indirect de-excitation rate, which is 
related to the total Lya scattering rate Pa (Field 1958) 



Pio=4P,/27, 



where 



cn^a^du, 



(3) 



(4) 



n^ is the number density of photons per unit frequency 
range, a{i') is the cross section for Lya scattering (Madau 
et al 1997). The brightness temperature in 21 cm is then 
(Field 1958, Chen & Miralda-Escude 2004) 



n = 25mK 
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It is worth mentioning here that recently Hirata & Sigurd- 
son (2007) showed that the spin temperature depends in 
general on atomic velocities Ts — Ts{v), if excitation and 
de-excitation coUisional processes are dominated by inter- 
atomic (H-H and H-He) collisions. This effect results in a 
sufBciently strong (up to 60%) increase of the 21-cm line 
width, however only in a 2% decrease of the emissivity. In 
this paper we focus mostly on the total emissivity, and there- 
fore neglect these effects. In principle, such a widening of the 



21 cm line can be important in the power spectrum of the 
brightness temperature on very small scales, X ^ 3 kpc. 



3 IONIZATION AND THERMAL HISTORY OF 
THE UNIVERSE 

Evolution of the fractional ionization is described by 

dXe 1 

"dT " {l + z)H{z) ' 



-{Rs{z) ~ Is{z) - Ie{z)) 



(6) 



where besides the standard recombination and ionization 
rates Rs{z), Is{z), the term Ie[z) is explicitly introduced due 
to presence of additional sources of UV photons. For the 
UHECRs this term can be written in the form (Peebles et 
al 2000, Doroshkevich & Naselsky 2002) 



Ie{z) = e{z)H{z)xii, 



(7) 



where t{z) = eo/(l+-z) is the production of ionizing photons, 
H[z) is the Hubble parameter, sh is the fraction of neutral 
hydrogen. For decaying particles this term is written as (CK) 



-fe(2) = X«/^rx- 



hvc 



(8) 



where Xi is the energy fraction deposited into ionization for 
which we made use the calculations by ShuU & van Steen- 
berg (1985), rup is the proton mass, fx = ^x{z)/Q,b{z), 
0.b{z) is the baryon density parameter, Q,x{z) is the frac- 
tional abundance of decaying particles, Tx is the decay rate, 
hvc is the energy of Ly-c photons. For short living particles 
Q,x{z) oc e~'"^' and their contribution to ionization and 
heating essentially vanishes at i > T^\ normally F^^^ is as- 
sumed to be comparable to the comoving Hubble time in 
the range of redshifts of interest. 

The gas temperature is determined by equation 



(1 + z)^ = 2r + kc „, , ,^ 2% ^ 

dz H{Z)[1 + fHe+Xc 



-{T-nmb) (9) 

K 



3kH{z){l + f„, + Xc) 

where the second term in the r.h.s. describes the energy ex- 
change between the gas and CMB photons kc = 4.91 x 
W~^^Tcmbj fue = 0.24, the third term is the heating from 
additional sources of the ionizing photons, K is the cor- 
responding heating rate which can be written in the form 
(CK) 



K = XhrripC ^Fj 



(10) 



where Xh is the energy fraction depositing into heating; as 
for Xi we used for Xh the results of ShuU & van Steenberg 
(1985). By order of magnitude x» ~ Xfe ~ 1/3 ^^r the con- 
ditions we are interested in. 

In the presence of UHECRs the final products are only 
Lya and Ly-c photons (Doroshkevich & Naselsky 2002). 
As pointed out by (Chen & Miralda-Escude 2004) the in- 
jected Lya photons change the gas temperature very little: 
Kc — 0. Ly-c photons have the energy slightly in excess 
of the hydrogen binding energy: ~ 13.6 eV-|-3fcr/2. There- 
fore the heating rate from the UHECRs produced Ly-c pho- 
tons Kc has to be calculated self-consistently together with 
equation (I10|l . However, as we will see later the gas kinetic 
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Figure 1. The kinetic (thin lines), spin (thick lines) and CMB 
(thick solid line) temperatures for the standard recombination 
e = (dashed), in the presence of UHECRs for the UV production 
rate e = 1 (dotted), and in the presence of decaying particles: long 
living with ^ = 3 X 10~^® s~^ (dot-dot dashed), and short living 
with Vx = 10"^5 s"\ fx{zeq) = 10"* (dash-dotted). 



temperature in presence of UHECRs is at most 150 K, so 
that the heating rate from Ly-c photons is always less than 
~ 2 X 10"^^e(l -I- zf^^ erg cm"^ s"\ The corresponding 
contribution to the gas kinetic temperature is less than a 
few percent. 

We apply a modified version of the code RECFAST 
(Seager et al 1999) to solve equations (|6]), pop . 



4 RESULTS 

Stellar and quasi-stellar sources of ionizing radiation be- 
gin to form at redshifts z < 20, and gas around them heated 
through photoionization can emit in 21 cm with a spot-like 
distribution. The spots have too small angular sizes, and 
can be detected only later, when at z ^ 15 star formation 
increases and forms pervading domains of sufficiently hot 
gas (Zaroubi & Silk 2004). Instead, UV photons from addi- 
tional ionizing sources, like UHECRs and decaying particles, 
illuminate and heat the IGM homogeneously, and therefore 
the signal from 21 cm can be more easily detected from the 
earliest redshifts. 

Fig. 1 shows kinetic and spin temperatures for a selected 
set of ionizing photon production rates by UHECRs and 
decaying particles. For the standard history (e = 0) the spin 
temperature is nearly equal to the CMB value at z ^ 20, and 
the intergalactic gas cannot produce distinguishible signal 
in 21 cm in this redshift range. The kinetic temperature 
for UHECRs with e = 1 is more than twice of the value 
in the standard ionization history (e = 0). The increase of 
kinetic temperature for higher e is due to an increase of the 




10-^ L- 



Figure 2. The contributions to the spin temperature from colli- 
sions (thick solid lines) and photon exitations (thin dashed lines) 
for the models with UHECRs upper panel for e = 0, 0.3, 1, 3 from 
bottom to top; long living particles mid panel: ^ = 6 X lO^'^^ s^^, 
3 X 10~26 s~i, 6 X 10-2^ s-\ 3 X 10-25 s'^ - from the low- 
ermost to the uppermost; i; = Xifx^x^ and short living lower 
panel decaying particles: Fx = 10"^'' s"'^, fxi^eq) = 0.5 X 10"^, 
Vx = 5 X 10-15 s-i, fxizeq) = 10-8, Fx = 10-15 s-1, 
fx(zeq) = 10-8, Fx = 10-15 s-1, fx(zeq) = 5 X lO"* - from 
the lowermost to the uppermost. 



fractional ionization of gas and a stronger coupling between 
the CMB photons and electrons. As a result, the gas kinetic 
temperature in this case shows similar variation with z as 
the CMB temperature. Due to coUisional de-excitations of 
the hyperfine structure level this increase in Tt unavoidably 
results in a decrease of the spin temperature, such that the 
difference between the spin and CMB temperature in the 
range z ~ 35 — 10 is for e — 1 greater than for e = 0; 
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Figure 3. The brightness temperature in the presence of UHE- 
CRs for the UV production rate e = (solid), e = 0.3 (dashed), 
e = 1 (dotted), e = 3 (dash-dotted). 



Ts — TcMB vanishes at z ^ 10. This means that the signal 
in 21 cm can be detected in the redshift range 2 ~ 10 — 35 
in absorption. The long living particles provide a permanent 
heating the injection energy rate 3kT/2 = K/{1 + fua + Xe), 
such that for a sufficiently high heating (ionization) rate 
shown in Fig. 1 the gas kinetic temperature in this case 
grows towards lower z as seen in Fig. 1 (thin dot-dot dashed 
line). Contrary, the short living particles inject heat with the 
rate T oc e~'"^', which manifests in a relatively fast decrease 
of the kinetic temperature at low z (thin dash-dotted line 
in Fig. 1). The HI spin temperature remains in both cases 
above the CMB temperature. 

As mentioned above, the heating of gas in models with 
UHECRs is small, and therefore the deviation of the bright- 
ness temperature from the standard value is mostly due to 
the Wouthuysen-Field (WF) effect from the UV photons 
produced by UHECRs. In Fig. 2 we show the contribu- 
tions of heating (collisional excitation and de-excitation) and 
Wouthuysen-Field effect to the spin temperature for UHE- 
CRs and decaying particles. 

Fig. 3 presents brightness temperature versus redshift 
for several values of e. An obvious qualitative difference be- 
tween the standard model (e = 0) and the models with 
e Js 0.3 is clearly seen: contrary to the standard case in 
all models Ti, flattens at 2 > 30 at the level T^ ~ -(10 - 15) 
mK. On the other hand at 2 < 25, where the standard model 
shows almost zero Tb all models with e Js 0.3 have bright- 
ness temperature between -5 and -10 mK. The signal of tens 
mK can though be easily swamped by a much stronger (of 
X^ 100 K) foreground emission at meter wavelengths. How- 
ever, multi-frequency observations seem to allow one to re- 
move the foreground, which is expected to be featureless in 
frequency space (Shaver et al. 1999, Di Matteo et al. 2002, 
Oh & Mack 2003, Gnedin & Shaver 2004, Zaldarragia et 



al. 2004). One can hope thus that an 1000 hours LOFAR 
(Low Frequency Array) and/or SKA (Square Kilometre Ar- 
ray) observation can discriminate between the standard and 
e ^ 0.3 models (see discussion in Sect. 4). 

Spatial variations of the brightness temperature con- 
nected with density perturbations in emitting gas are to 
be more easily distinguished in the redshifted 21 cm emis- 
sion. For this purpose we calculated the derivative of the 
brightness temperature Tb over the amplitude of the den- 
sity fluctuation S, where the perturbed density field is as- 
sumed in the form p — po{z){l + S), with po{z) being the 
background baryon density. We calculated the deviations of 
the spin and kinetic temperatures starting from zq — 1000 
and followed the evolution of perturbed regions according 
to equation (|10[1 to the redshift of interest z, accounting 
the dependence of kinetic coefficients fcio and je on tem- 
perature. In Fig. 4 we show results for the standard recom- 
bination model: the upper panel of Fig. 4a illustrates the 
deviation STb of the local brightness temperature from the 
all-sky one for a separate perturbation of a given density 
amplitude 5. As expected, the overdense regions show larger 
amplitude in absorption. The dependence of STb on redshift 
(x-axis) weakens at higher z. For the sake of clearity in the 
lower panel of Fig. 4a we show also the derivative dTb/dS 
versus redshift and amplitude: at z = 20 — 40 it weakly de- 
pends on the amplitude S, although changes very sharp ver- 
sus redshift. The total signal comprises contributions from 
the whole spectrum of perturbations and is represented by 
the temperature power spectrum (Barkana & Loeb 2005ab, 
Hirata & Sigurdson 2007) 

Pt, = P^o{k)+ii^P^2{k)+p^P^4k), (11) 

where ^ = fc||/fc is the cosine of the angle between the line 
of sight and the wavevector. 



f dTb 
\ dS 



= kTb-^^Ps 
do 



(k), 



(12) 



.(fc), 



Tb is the mean brightness temperature, Ps{k), Pv(k) axe 
the power spectrum of the baryon density and velocity 
fluctuations, Pi,^(fc), their cross-spectrum. Fig. 4b shows 
5Tb{k,z) = wWPii/2tv^ for the standard model in the 
plane {z,k): the dominant contribution stems from the 
density fluctuations, the contribution from density-velocity 
cross-spectrum is always around 20-30 %, the velocity fluctu- 
ations contribute less than 10-20 % at 2 ^ 40, and increases 
(up to 30%) at 2 ~ 50; these estimates are given for the 
average (over 47r) value of p'^ = 1/3. 

In Fig. 5 we show the distributions STb, dTt/dS and 
5Tb{k, 2) for the model with UHECRs e = 1. As in the stan- 
dard model 21 cm line can be observed in absorption, how- 
ever there is a clearly seen distinct feature between the two 
cases: the standard model shows gradient of the temperature 
power spectrum from low amplitudes and redshifts (right 
bottom corner in Fig. 4b) towards higher amplitudes and 
redshifts (left upper corner), while the model with UHECRs 
reveals an opposite behaviour of the gradient of STb{k, z) - 
from the left lower to the right upper corner, and in addi- 
tion, a much smaller gradient in redshift (i.e. in frequencies) 
at lower amplitudes. 
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Figure 4. a) The diflference between the local and all-sky bright- 
ness temperatures (ST;, = Ti,(5) — Tj, (upper panel), the derivative 
dTf,/dS (lower panel) as functions of redshift and density variation 
5; b) the temperature power spectrum STi,{k, z) as a function of 
redshift and wavenumber k: standard model. Numbers on isocon- 
tours are given in mK. 



Figure 5. Same as in Fig. 4 for model with UHECRs: e = 1. 
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Figure 6. Solid lines — the brightness temperature in the pres- 
ence of long Hving particles Vx <C Hit: ^ = 6 X lO^^'^ s^^, 
SxlO"^'^ s"\ exlO"^*' s^l, 3x10^25 s^i -from the lowermost 
to the uppermost; ^ = Xifx^x- Dashed lines - the brightness 
temperature in the presence of short living particles Tx X^ ^o: 
Vx = 10"" s-i, fx(z^q) = 0.5 X 10"*, Tx = 5 X W^^ s''^, 
fx{ze<i) = 10-8, Fx = 10-15 s-1, /x{ze,) = IQ-*, Tx = 
10~15 s~l, fx{zeq) = 5 X IQ-* - from the lowermost to the 
uppermost; the thick part of the solid curve and the dot-dashed 
thick parts of the dashed curves indicate the predominance of WF 
effect. 



The models with decaying dark matter can produce 
21 cm in emission and in absorption as well. Fig. 6 shows 
the all-sky brightness temperature for the long living (solid 
lines) and short living (dashed lines) decaying particles. The 
gas kinetic temperature grows with heating rate and can be 
lower or greater than Tcmb depending on ^. Therefore, at 
low heating rate (^ ^ 6 x IQ-^^ s-^) gas can be observed in 
21 cm only in absorption with a lower absolute brightness 
temperatue than in the standard case, essentially approach- 
ing it due to decrease in the heating rate (compare the lower 
solid line in Fig. 6 with the solid line in Fig. 3). The tran- 
sition from absorption to emission occurs at ^ Js IQ-^® s-^. 
The long living decaying particles affect the spin and the 
brightness temperature mostly through collisional heating 
(as seen in Fig. 2) , however in some cases Wouthuysen-Field 
effect can become more important: this is the model with 
^ = 3 X lO"^'^ s"^ at z > 40 - the contribution from WF 
effect is shown in Fig. 6 as a thick part of the corresponding 
curve. 

Brightness temperature for short living particles Tb in- 
creases with the heating rate similar to what occurs for the 
long living particles, and correspondingly at a sufficiently 
high heating 21 cm can be observed in emission. There is, 
however, a qualitative difference between the dependence 
of brightness temperature on redshift in these two cases: 
for the long living particles Tb{z) has always negative sec- 



ond derivative, while in the case of the short living particles 
Tb{z) has an inflection point. Note in this connection that 
Tt{z) curves in models with the UHECRs have always pos- 
itive second derivative. Similarly to the models with long 
living decaying particles, for short living particles a mostly 
collisional (heating) contribution to Ts and Ti, is reversed to 
predominantly Wouthuysen-Field effect for the models with 
Tx = 10-15 s-i and fx{zeq) = 10"^ and with Fx = 10"^" 
s-i, fx{zeq) = 5 X IQ-* at z > 30 and z > 20, respectively; 
this is also shown by the thick parts of the corresponding 
curves in Fig. 6. The inflection point on Tb{z) nearly corre- 
sponds to redshifts where WF effect becomes dominating. 

The differences between the second derivatives of the 
all-sky brightness temperature Tb{z) is reflected to a certain 
extent on the behaviour of the temperature power spec- 
tra on (z, k) plane. As seen in Fig. 7b the lines of equal 
spectral power for the long living particles reveal a smooth, 
nearly uniform, gradient over the whole range of redshifts 
and wavenumbers (amplitudes). At the same time, for the 
short living particles (Fig. 8b) these lines thicken at the low 
z side, with a stronger gradient, similarly to the model with 
UHECRs where d'^Ti,/dz'^ is always positive. The behaviour 
of the temperature difference 5Tb and its derivative dTb/dS 
on the (z, k) plane are also distinct for these two models of 
decaying particles (Fig. 7a and 8a): the short living parti- 
cles show thickening near the low end of redshifts. These 
differences though are not well pronounced (as, for instance, 
between the standard model and the model with UHECRs 
shown in Fig. 4 and 5), and it seems to be challenging to 
reveal them. In the high redshift end z > 20 (lower fre- 
quences), the STt{k,z) maps for the two models look quite 
similar. However, at z < 20 (higher freqiences) the tem- 
perature power spectrum at a given wavenumber decreases 
towards lower redshift faster in the model with short living 
particles, which can be recognizable in the frequency space. 

The differences between the second derivatives of the 
all-sky brightness temperature over redshift d^Tb/dz^ for the 
models with UHECRs and long- and short living particles, 
or the corresponding differences between the temperature 
power spectra, imprint in spectral features of the 21 cm 
line. This circumstance may have a principal significance 
for choosing a strategy for observational discrimination be- 
tween manifestations from the three sources of the ioniz- 
ing photons, provided that observations are possible for a 
wide range of redshifts, from z = 10 to z = 50. However, if 
the qualitative difference between Tb{z) and 5Ti,{k, z) curves 
for the UHECRs and long living particles seems to be rela- 
tively easily observationally distinguishible, it looks prob- 
lematic for the curves Ti,{z) and STb{k,z) in cases with 
the long and short living particles because the inflection 
of Ti,{z) in the latter case lies within AT[,(z) ~ 5 mK, 
and the differences between the 5Tb{k,z) maps are basi- 
cally within a few mK for wavenumbers from, e.g., 0.1 to 
1 Mpc-i. From this point of view "two-color" diagrams, 
connecting the differences between the brightness temper- 
atures at different wavelengths, can be complementary to 
the analysis of Tb{z) curves and ST^ik, z) maps. We show an 
example of such "two-color" diagrams for the three sources 
of UV photons: UHECRs, long and short living decaying 
particles in Fig. 9. Specifically, we plot the relative all-sky 
temperature differences A^T23 — AT(,(A2) — Arj,(A3) versus 
A^Ti2 = ATb{\i) - ATb{X2), where Ai = 210 cm, A2 = 420 
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Figure 7. Same as in Fig. 4 for the model witli long living par- 
ticles Tx < -ffo for 5 = 6 X 10-26 s-i. 



Figure 8. Same as in Fig. 4 for the model with short living 
particles Vx ^ Hq: Fx = IQ-i^ g-i, fxizeq) = IQ-*. 



Yu. A. Shchekinov and E. 0. Vasiliev 



CM 



15 


1 1 


1 1 


1 1 




1 1 


10 


- 








♦ • ■ 


b 


- 








•- 


n 


- 






E 


' • 


-5 


- 




D 




♦ : 




- 




D 




. 


1U 


- 








- 


15 


- 








- 


20 


- D 








- 


25 


- 








- 


30 


- 


o 






- 


35 


- o 

- 1 1 1 


1 1 


1 1 




1 1 



-16 -12 

.2. 



-8-4 4 



Figure 9. "Two-color" diagram for the models with UHECRs 
and decaying particles. Open symbols show differences in bright- 
ness temperatures in emission, filled - in absorption; half-filled 
correspond to emission at low z (lower half is open) and absorp- 
tion at high z (upper half is filled), and vise versa. Circles show 
(A'^ri2, A^T23) for UHECRs, diamonds - for long living decaying 
particles, squares — for short living particles, the triangle depicts 
the standard (no additional ionizing photons) case. 



cm, A3 = 840 cm, for a set of models for these three cases, 
here /S.Ti,(\i) is the difference between the total signal at 
\i = 21(1 -I- Zi) cm and the temperature of the foreground 
emission at this wavelength. Since the foreground emission 
has temperature of several orders of magnitude higher than 
the expected from 21 cm, determinantion of the 21 cm signal 
with using such a substraction procedure seems challenging. 
However, as pointed out by Di Matteo et al. (2002), and Oh 
& Mack (2003), the prospectives for measurements of the 
all-sky signal from redshifted 21 cm are not totally dim be- 
cause the foreground emission is expected to be featureless 
in frequency (see also discussion in Zaldarriaga et al. 2004, 
Sethi 2005). The three cases fall into quite distinct regions 
in the "two-color" plane with a separation of A'^T > 5 mK, 
and thus can in principle be discriminated observationally. 
The parameters of decaying particles and UHECRs we 
explored in this work are consistent with the current con- 
straints. CK found using WMAP data ^ < lO"^'* s"^ for 
the long living particles, and a slightly weaker upper limit 
for the short living ones. Recent results from WMAP give a 
lower value of the optical depth r ^ 0.09 (Spergel et al 2006). 
This strengthens the constrains on the parameters of decay- 
ing particles compared to those inferred from the WMAP 
data of the first year. However, as pointed out by CK de- 
caying particles with short lifetimes do not affect singnifi- 
cantly the optical depth and are thus less constrained. On 
the other hand, the upper limit for ^ = 3 x 10^^^ s^^ for 
long living particles, we consider here, lies within the re- 



strictions on possible ionization rates corresponding to the 
optical depth obtained from the third year WMAP data. 
For UHECRs Doroshkevich & Naselsky (2002) inferred from 
the MAXIMA-1 and BOOMERANG data (de Bernardis et 
al 2000, Hanany et al 2000) e < 3. Observations of HI at 
z = 10 — 50 in 21 cm can provide further constraints. 

The sensitivity of ongoing long- wavelength experiments, 
such as LOFAR, 21CMA (former PAST), MWA, LWA, SKA 
and LUDAR (Carilli 2006) seems to be sufficient to detect 
signal in 21 cm affected by decaying particles and UHECRs 
at pre-ionization epochs. The minimum background source 
flux density reqiured to detect an absorption feature with 
the signal to noise ratio S/N is (Furlanetto 2006) 



Omin — tb - I 

V 5 ; 



10 



-3 



lO^m' 



T,, 



400 K 



^ J \ '^ J \ ^■•'ff 
100 kHz 1 weck \"-^ 
Aiy ti / ™ ^ 



(13) 



where r is the optical depth in 21 cm, Aeff is the effec- 
tive area of telescope, Az^ is the bandwidth of channel, ti 
is the total integration time; in the low frequency (~ 100 
MHz) range the system temperature Tsys is dominated by 
the Galactic synchrotron background and scales as z^~^'^, or 



in terms of redshift (Chen 
etto et al. 2006b) 



2000 



21 



K. 



Miralda-Escude 2006, Furlan- 



(14) 



For a conservative estimate we take the bandwidth of a chan- 
nel Av = 4 MHz, as a typical value for future telescopes. 
The effective area of the telescopes varies from ~ 3 x lO'^m^ 
for LOFAR (at 75 MHz) to 10*^ m^ for SKA and LWA (for 
LUDAR it can reach up to ~ 10^ m^). The optical depth 
in 21 cm increase from 10"'^ at z ~ 10 to 0.01 at 2 ^ 15 
(Furlanetto 2006). Therefore, for the redshifts 2; ^ 10 we are 
interested in r can be taken > 10^^. The integration time is 
assumed to be of 10 weeks. Thus, the minimum background 
flux is 12 mjy at z = 20 - 40 for LOFAR, 0.4 mjy at z = 20 
for SKA and LWA, and 2.3 mJy at 2 = 40 for LWA. In the 
case of UHECRs the difference of a smooth signal expected 
between the models with e > 0.3 — 3 is > 5 mJy, and is even 
larger in the models with decaying dark matter. This sup- 
ports the conclusion that observations on future telescopes 
can discriminate between the standard model and those af- 
fected by UHECRs and decaying particles. 



5 CONCLUSIONS 

In this paper we have considered the influence of UV photons 
from decaying dark matter particles and UHECRs on the 
ability of neutral hydrogen to emit and/or absorb in 21 cm 
at redshifts z = 10 — 50. We have found that 

• the three sources of additional ionizing photons: long 
living and short living unstable dark matter particles and 
UHECRs produce fairly distinct dependences of the all-sky 
brightness temperature on redshift Tb{z) - the flrst and the 
third give negative and positive second derivatives of the 
curves Tb{z), while the second has Tt{z) with an inflection 
point. Although the all-sky 21 cm emission is expected to 
be swamped by the owerwhelming foreground signal, one 
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can hope that due to the lack of features in frequency space 
the latter can be removed. Moreover, the maps of the power 
spectrum of the brightness temperature on the wavenumber- 
redshift plane, reveal clear differences (seen in gradients of 
the temperature power spectrum over redshift) between var- 
ious models of the particles. 

• these features manifest in the frequency space of 21 
cm line. From this point of view three wave-band obser- 
vations at Ai, A2 and A3 and a "two-color" diagram for 
the relative (with respect to the foreground emission at a 
given wavelength) all-sky temperature differences A^T23 = 
ATi,(A2) - Ar6(A3) versus A^Tiz = Ar6(Ai) - AT6(A2) can 
provide an additional tool for discrimination between the 
sources of ionizing photons in the end of dark ages. 

In general, decaying dark matter particles can have 
strong effects on overall history of the universe: they 
may change its thermal evolution and cosmological nucle- 
osynthesis (Scherrer 1984, Vayner et al 1985, Vayner & 
Shchekinov 1986), dynamics of large scale structure forma- 
tion (Doroshkevich et al 1989, Bharadwaj & Sethi 1998, 
Cen 2001), reionization regime (Sciama 1982, Dodelson & 
Jubas 1994, Hansen & Haiman 2004, Chen & Kamionkowski 
2004, Kasuya et al 2004, Kasuya & Kawasaki 2004, Pier- 
paoli 2004, Mapelli et al 2006), and even formation of the 
fisrt stars though enhancement of H2 molecule formation 
(Shchekinov & Vasiliev 2004, Vasiliev & Shchekinov 2006, 
Biermann & Kusenko 2006, Ripamonti et al 2007, Furlan- 
etto et al. 2006a). 21 cm emission from the epochs ending 
the dark ages can carry the imprints from decaying particles, 
and seems a promising tool for understanding their prop- 
erties. Future radio telescopes (such as LOFAR, 21CMA, 
MWA, LWA and SKA) seem to have sufficient flux sensitiv- 
ity for detection the signal in 21 cm influenced by decaying 
particles and UHECRs. 
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